Absence of two energy scales in the two-impurity Kondo Model 
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It is believed that the successive antiferromagnetic scattering of the conduction electrons on two 
magnetic impurities in a metal, induces a magnetic interaction between the impurities which sets 
an energy scale for the system, in addition to the Kondo temperature. However, it is shown here 
that this contradicts a scaling law for the two- impurity Anderson model in the magnetic limit which 
becomes exact for a large bandwidth. 



In heavy fermion materials localized, strongly correlated / electrons interact with delocalized, weakly correlated 
conduction electrons. Experiments show that above a temperature T* the magnetic moments of the / electrons behave 
asymptotically as free. 

Well below T* they can either be pheiiomenologically described as a Fermi liquid, or magnetic correlations between 
, the local / moments become dominant .El 

When searching for a microscopic mechanism, one is guided, by the behavior of dilute magnetic alloys which are 
well described by the single-impurity Anderson HamiltonianB At low temperatures the magnetic moment of the 
impurity is screened by particle-hole excitations of the conduction band. The system behaves universally, that is, all 
observables scale with one energy UbTk where Tk is called the Kondo temperature. On the other hand, the heavy 
fermion materials have, in orders of magnitude, one magnetic "impurity" atom per unit cell. The hypothesis for such 
I ' systems is that it can be described by the Anderson-lattice Hamiltonian. 

It is far from obvious why the single-impurity Kondo effect should survive in the lattice model. It is believed 
that magnetic interactions between the local / moments compete with the Kondo effect. This magnetic interaction 
between the local / moments is supposed to be induced by successive scattering of the conduction electrons on the two 
i^i - impurities, and called Ruderman-Kittel or RKKY interactionJjTJ Doniacbo estimated how much energy the system 
gains if it orders antiferromagnetically. By second-order perturbation theory in the exchange coupling J between the 
1 ■ conduction and the / electrons this RKKY energy was found to be of the order 
> ■ 
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where p is the density of states at the Fermi energy of the conduction band of width D. Doniach compared it with 
the exponentially small energy gain 
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if the magnetic atoms are considered to be independent magnetic impurities. Depending on the size of J, it could then 
I be estimated whether the material is in a magnetic or nonmagnetic phase, at low temperatures. This simple picture 
I ■ deals essentially with a two-impurity Hamiltonian, the standard version of which is the two-impurity Anderson or 
Kondo model. It has had great impact on the experimentalist's point of view. In fact, standard reviews quote the 
Doniach picture to explain, at least qualitatively, the respective experimental observations.tm 

In this article it will be shown that this simple picture cannot be valid. To this end let us discuss firstly the various 
approaches to the two-impurity problem. _ _ 

For the two-impurity Kondo model, summing of leading-order divergent diagramaD in the spirit of Abrikosovcl 
showed that at decreasing temperature, the RKKY interaction as well as the Kondo interaction are renormalized, and 



5_j ■ eventually diverge at low temperatures. Hence the picture drawn by Doniach would be invalid. 

This can be contrasted with results from poor man's scaling:Bt] Depending on the distance between the impurities, 
the RKKY interaction / may either be ferro- or antiferromagnetic. If / is ferromagnetic and larger than ksTx, it is 
predicted that the magnetic moments of the two impurities first form an effective spin one and are then quenched 
in a two-stage process, as temperature decreases. On the other hand, a large antiferromagnetic RKKY interaction 
suppresses the Kondo effect and the two impurity spins form a singlet for k B T < I. It is then argued that / sets an 
energy scale for the system at which the RKKY interaction could be observed. 

This analysis was-s«eported by a more recent calculation with the help of a variant of the numerical renormalization 
group by Silva et aLl!30 (in which the Nth hopping matrix element of the half-infinite chain to which the two impurities 
couple is vanishing as A~ N with A = fO, and subsequently the results of different band-widths D/A e , e = . . . 1 are 
averaged over). In fact, those results imply that the Doniach picture is valid. 
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However, those results contradict others£3 using the conventional renormalizati on g roup scheme (A = 3) where only 
one low-energy scale was seen. What is more, the approaches of Refs. and [l(]are inconsistent. To clarify the 
point, only the case of an antiferromagnetic RKKY interaction is discussed. Silva et al. interpreted their data with 
the help of poor man's scaling. They saw that the magnetic moment of the impurities is quenched as temperature 
decreases below I/ks, at which temperature / and J themselves are not much renormalized. They claimed that 
therefore the two impurities would form a singlet for fc^T < /, consequently the renormalization of J would stop, and 
the low-temperature behavior could be determined from perturbation theory, the lowest nonvanishing order of which 
is oc J 2 /I. A simple estimation shows however that this perturbative expansion is not allowed: Because p oc 1/D, 
the RKKY interaction IJs always a factor J/D smaller then the Kondo interaction J. This can be read off their 
numerical results as wellJiJ 

Incidentally, if the arguments given by Silva et al. were correct they could be repeated for the antiferromagnetic 
single-impurity Kondo model: At a temperature of the order of the coupling constant J the impurity spin and an 
electron of the conduction band would form a singlet, the renormalization would stop, and one would be left with a 
"Kondo effect" at an energy scale of the order of J, which is certainly not true. 

To summarize, the Doniach picture and its subsequent refinements as described above, imply that the renormaliza- 
tion of J stops as temperature becomes lower than the RKKY coupling. This RKKY interaction then sets an energy 
scale for the system's response, at least for small coupling J <C D. Therefore it would be helpful if the exact energy 
scales of a metal with two magnetic impurities could be calculated in the universal limit of large D. The method 
which is usually employed for that purpose is the diagrammatic renormalization group. The major obstacle is the 
solution of the Dyson equation connecting the dressed propagators and their self-energies which can be done only 
perturbatively. However—the low-energy behavior is beyond perturbation theory. An explicit solution of the Dyson 
equation can be avoidedEj by means of the variational principle of Luttinger and WardEJ~t3 with the help of which 
the energy scales will follow. 

The two-impurity Anderson model serves here as the standard model for two magnetic impurities in a metal. 
Mean-field|-approaches for this model have the drawback that the RKKY interaction appears only as a higher-order 
correction, E3 and therefore a magnetic interaction between the impurities had to be added by hand, giving a model 
with an adjustable strength of the magnetic exchange between the impurities. 

Here, however, a direct exchange interaction between the impurities will not be taken into account, because the two 
impurities model two localized / orbitals on different lattice sites, with virtually no overlap, and it is the purpose of 
this article to investigate the role of their induced magnetic interaction, that is, the RKKY interaction. 

The two impurities are assumed to sit at (0,0, ±i?/2) symmetrically to the origin of the coordinate system, along 
the z axis. The Hamiltonian H is then invariant under the parity transformation. Hence the creation operator for a 
conduction electron of momentum p and | smu m and for an impurity electron can be expressed in terms of their even 
(a = 1) and odd (a — —1) combinations ,[3Ej 

dt — — — ( ct + (JC + - \ 
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Here ct m creates a conduction electron with internal quantum number m = 1 . . . N, momentum p, and energy e(p). 

If |vac) is the vacuum, then 6 + |vac), /^|vac), and d+„|vac) denote the unoccupied, singly occupied, and doubly 
occupied impurity configurations, with % = 1,2 labeling the impurities. The energy difference between the singly 
occupied and unoccupied as well as between the doubly and singly occupied impurity configuration is e/. ff m is a 
fermionic operator, and b + and d^ nn are bosonic. Double occupancy at each impurity site is suppressed by requiring 
b+b + Y.am famfam + E- ™>„ d trnn d arnn = 1 - The unoccupied configuration 6+|vac) has even parity. The effective 
Hamiltonian then reads,t£l 

H = H C + H f + H 1 , 

^ ^ C emcr C emcn 
ea 



m>n,cr 
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+ 77f/V E r ^(e)(4T, m n/ ra C Tem + ff ' c -) 

+ :7TF E ^C 6 ) ( d -l, mm f-*mCcT,e m + H.c) . (2) 

The impurity hybridization with the conduction band is proportional to V. The Boltzmann constant is set to unity 
so that temperature is measured in units of energy. The dispersion can be assumed as isotropic. The density of states 
p is assumed to be constant, and the such linearized dispersion is cut off at ±D. T-h|C effective hybridization matrix 
element can then be expressed in terms of the conduction band density of states as 



VAe) = VVPB(1 - ^ 1 + a ^ - ^ g) , (3) 

with Icf the Fermi wave number and O the step function. Thus the system has been reduced to an effective one- 
impurity Hamiltonian, the additional parity quantum number keeping track of the original two impurities. 

The one-particle impurity propagators Rf for the unoccupied (/ = 0), singly occupied (/ = am), and doubly 
occupied (/ = amn) impurity configurations, determine the impurity part Zf of the partition function via a line 
integral, the path of integration encircling all poles of the integrand 

where TV/ denotes the tp.ce over the impurity configurations. The propagators Rf can be calculated by the standard 
diagrammatic technique^ which isrfeviewed in Ref . 20 . The vertices and naked propagators are shown in Fig. |l|. 



Within the variational principle,!^ a functional T of the propagators Rf is defined in terms of skeleton diagrams. 
For the impurity part of the partition function, T is given by 

T = ^ f Ea - £)4 B) (*/(*))*/(*) 

J n 

+Hz-H f -j2z f n \Rf(m}- ( 5 ) 

n 

Here denotes the nth order self-energy of Rf, expressed in terms of skeleton diagrams. At the saddle point 
with respect to variations of Rf, the functional T equals Zf and the Dyson equation holds as a self-consistencjz. 
equation.li3 So T depends on parameters such as e/ explicitly only via Hf, and not implicitly via the propagators .113 
In addition, T depends explicitly on V only via the prefactor V 2n of the 2nth order self-energy £^ 2 ™\ To determine 
the T dependence, the internal integration variables z and e, as in Eqs. (|) and (@) are replaced by Tz, Te, and D/T, 
respectively. The variational principle remains unaffected. T depends now explicitly on T via the prefactor T n of 
a skeleton diagram of 2nth order, the term Tz in the logarithm, and D/T because of Eq. (|^). Similar to Ref. |l3| it 
follows for the impurity part Ff of the free energy 

Approximations fulfilling the variational principle can be generated by using a subclass, the so-called families of 
skeleton diagrams.Eil If, for instance, only the family of second-order skeleton diagrams is kept in T, the self-energies 
are (/ denotes the Fermi function), 



E^UnC*) = ^ E / def(e)V T 2 (-e)[R an (z + e) 

err— const 

+R am (z + e)], 
4t(*) = ^ fdef{t)[\ V r 2 (e)R„ T , mn (z + e) 
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+V 2 (-e)R Q (z + e) + V 2 a (e)R- hmm (z + e)] , 

a.ra 

Universal behavior of the system manifests itself when D becomes larger than all other energy scales of the system. 
T depends explicitly on D only via V 2 as in Eq. (||), 

Dd D T = -Tr f I ^e-t*Dd D E ^f\z)R f (z). (8) 

n 

The skeleton diagrams for YrpRf are shown in Fig. ||. 

The effective density of states of Eq. (||) scales with D as a function of e and is finite at the Fermi energy e = 0. 

Hencejiij as in the case for a density of states with a sharp cutoff at ±D, differentiating a skeleton diagram for S^ n) i?/ 
with respect to D amounts to removing one curved conduction electron line and replacing the internal propagator by 
its value at the cutoff oc 1/D. Therefore a skeleton diagram of second order contributes oc D/D to the logarithmic 
derivative. 

A skeleton diagram of higher than second order contains vertex corrections; hence there lie under each conduction 
electron line at least two impurity propagators, because otherwise this diagram would have a self-energy insertion and 
noJ— be a skeleton. Its contribution to the logarithmic derivative is therefore cx D/D 2 and can be neglected for large 

DB 

Therefore, Eq. (||) only the skeleton diagrams for S^iZ/ are necessary to obtain the exact energy scales of the 
system in the universal limit of large D. A spectral decomposition of the impurity propagators yields, 



°dD Ff = ~l ^Z~\ 2 V ° (0) [/w(w) + *>W 

crTmn 



Together with the identities 
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Z f —F } = e- flu du$2p <rm (u) + 2 E fcmmM], 
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the following scaling equation is obtained, 



D d3 Ff = py2 (1 " 1/N) fe~ f Ff ~ 2py2 ' (9) 

Eqs. (0) and m) imply the scaling law 

F f -E = Tg(T, V, e f , D, p) = Tg (J, (10) 

in the magnetic limit — e/ ^> pV 2 , with T = irpV 2 /N. In this limit, the Anderson model is equivalent to the two- 
impurity Kondo modelEl at low temperatures T <C pV 2 , with Tk = D pV 2 / D exp [ey /pV 2 ~\ . By the same meansEl 
it can be shown that in this limit all observables scale as in Eq. (px|). 

Hence the two-impurity Kondo model has in the limit of large D only one low energy scale, which is proportional 
to the single-impurity Kondo temperature. Observe, however, that the proportionality factor cannot be fixed within 
this approach and will of course depend on kpR. 
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If the RKKY interaction would set additional energy scales as discussed in the introduction, then the scaling law 
would have the form, 



f T T T T\ 
F f -E = g (11) 

\ll< I J-KA J- + J 

where Tka an d T+ are the temperatures of the two-stage Kondo effect for the ferromagnetic RKKY interaction]!^ 



This, however, contradicts the scaling law (10) above. 

To conclude, if the high-energy cutoff of the metal is so large that terms of the order T/D or J/D can be neglected, 
there is for two magnetic impurities in a metal only one low-energy scale. The induced indirect exchange interaction 
between the local / moments which is mediated by successive scattering of conduction electrons, will be strongly 
renormalized at low temperatures. However, the scaling equations alone do not predict the nature of the ground state 
for two magnetic impurities in a metal. In view of Eq. mq) it seems doubtful whether a competition between the 
formation of magnetic order and a Fermi liquid of heavy quasiparticles can be understood without adding explicitly 
a magnetic interaction between the magnetic / moments. 

While preparing this work the author has benefited from numerous discussions with Tom Schork, Karen Hallberg, 
Karlo Penc, and Professor Peter Fulde. 
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FIG. 1. Vertices for the two-impurity Anderson model. A dashed, singly or doubly arrowed line represents the naked 
propagator of the singly or doubly occupied impurity configuration. A single a or r labels the parity ±1, and err their 
product, m or n labels the magnetic quantum number, and mn a pair of magnetic quantum numbers for the double occupied 
configuration. A wavy respective solid line represents the naked propagator of the unoccupied impurity configuration or the 
conduction electron. The product of parities of incoming particles equals that of the outgoing ones. 



FIG. 2. Skeleton diagrams of second order for the two-impurity Anderson model 
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